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1.  Introduction 

The  conservation  laws  of  isothermal,  isentropic  or  adiabatic 
thermoelasticity,  in  all  their  standard  variants  (Lagrangean  or 
Euclidean  formulation,  solids  or  fluids,  one  or  several  space 
dimensions,  etc.)>  lead  to  systems  of  quasilinear  hyperbolic  equations,  j 
A  feature  of  such  systems  is  that  the  Cauchy  problem  does  not  have  global 
smooth  solutions,  even  when  the  initial  data  are  very  smooth,  due  to 
the  formation  of  shock  waves.  However,  global  solutions  exist  in 

the  class  of  functions  of  bounded  variation* /in  the  sense  of  Tonelli- 

/  ; 

Cesari  [8]. 

When  the  material  is  viscous  and/or  heat  may  diffuse,  dissipative 
mechanisms  emerge  in  thejsystem  of  conservation  laws,  which  manifest 

i 

themselves  through  the  appearanceof  "parabolic”  terms  or  "memory" 
terms.  The  same  phenomenon  also  arises  in  the  context  of  the  theory 

I 

of  /chemically  reacting  media  with  dissipation  induced  by  diffusion. 

^  A  dissipative  mechanism  may  affect,  in  general,  the  asymptotic  j 

behavior  as  well  as  the  smoothness  of  solutions.  Ranked  according  to  j 

\ 

effectiveness,  dissipative  mechanisms  may  be  classifed  into  ftrj  those  j 
which  are  so  powerful  as  to  smoothen  out  even  rough  initial  data,  always 
yielding  smooth  solutions;  those  that  preserve  the  smoothness  of 
smooth  initial  data  but  are  incapable  of  smoothening  rough  initial 

t  i 

data;  (c)  those  that  preserve  the  smoothness  of  smooth  and  "small" 

initial  data  but  cannot  prevent  the  breaking  of  smooth  waves  of 

, ,  $  : 

large  amplitude;  (d)  those  that  are  not  capable  to  prevent  even  the 

breaking  of  smooth  waves  of  small  amplitude. 

From  the  viewpoint  of  analysis  as  well  as  continuum  physics 

it  would  be  useful  to  classify  the  conservation  laws  for  the 


standard  material  classes  into  one  of  the  above  categories.  Although 
a  number  of  scattered  results  in  this  direction  are  already  available , 
there  are  still  many  unanswered  questions  and,  more  importantly,  there 
has  been  no  attempt  to  place  the  existing  information  into  a  general 
framework. 

In  this  lecture,  I  shall  outline  a  research  program  towards 
understanding  the  role  of  dissipation  and  estimating  its  effectiveness. 
I  will  exhibit  a  (far  from  exhaustive)  list  of  representative  examples 
and  I  will  discuss  what  has  been  established  and  what  has  been  con¬ 
jectured,  for  each  case. 

2.  Complete  Parabolic  Damping 

The  simplest  example  of  an  equation  with  complete  parabolic 
damping  is  provided  by 

(2.1)  ut  ♦  f(»)x  -  uxx, 

(f(u)  nonlinear  and  smooth)  whose  behavior  is  to  be  contrasted  to 
that  of  the  hyperbolic  conservation  law 


(2.2)  ut  +  f^x  *  °* 

The  Cauchy  problem  for  (2.2),  with  initial  data  u(x,0),  of  bounded 
variation,  admits  a  solution  in  the  class  BV  of  functions  of 
bounded  variation,  in  the  sense  of  Tonelli-Cesari .  No  gain  would  be 
made  by  assuming  that  u(x,0)  is  smoother,  even  analytic!  In 
contrast,  the  initial-value  problem  for  (2.1)  with  rough,  just 


bounded  measurable,  initial  data  always  has  a  smooth  solution. 

This  type  of  behavior  is, of  course,  due  to  the  power  of  parabolic 
damping  and  characterizes  a  broad  class  of  systems  of  equations  in 
the  form 


(2.3)  wt  +  f(w)x  «  Awxx 

where  A  is  a  positive  definite  matrix  and  f  satisfies  certain 
technical  smoothness  and  growth  assumptions. 

3.  Incomplete  Parabolic  Damping 
In  this  section  we  will  discuss  the  representative  systems 

(3.1) 

(3.2) 

C3.3) 


(3.1)  are  the  equations  of  motion  of  one 'dimensional  viscoelastic 
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materials  of  the  rate  type,  while  (3.2)  and  (3.3)  express  the  con¬ 
servation  laws  of  momentum  and  energy  of  one-dimensional  thermo¬ 
viscoelastic  and  thermoelastic  materials,  respectively. 

The  physically  natural  assumptions  are 


(3.4)  Pu(u,8)  <  0,  e0(u,0)  >  0,  eu(u,6)  =  9  2  [P-lufr9J  ] 

(3.1)  should  be  contrasted  to  the  system  of  equations  of  motion 
of  one-dimensional  elastic  materials, 


(3.5) 


ut  -  vx  -  0 


vt  +  p(u)x  “  °» 


while  (3.2)  and  (3.3)  should  be  contrasted  to  the  conservation  laws 
for  adiabatic  processes  in  one -dimensional  thermoelasticity: 


(3.6) 


Ivt  +  p(u,8)x  *  0 

[e (u,0 )  ♦  ^]t  +  [p(u,0)v]x  *  0. 


(3.5)  and  (3.6)  are  typical  examples  of  systems  of  hyperbolic 
conservation  laws. 

System  (3.1)  can  be  written  in  the  form  (1.3)  where,  however, 
A  is  only  positive  semidefinite  (whence  the  term  "incomplete 
parabolic  damping").  Intensive  investigation  for  System  (3.1) 
(e.g.  [1,3,9,10])  has  generated  a  chain  of  existence  theorems 
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revealing  several  function  classes  whose  smoothness  is  preserved 

00  2 

by  solutions.  Typical  examples  of  such  classes  are  L  *  L 
(cf.  [1])  and  C^+°  *  C^+a  (cf.  [3]).  However,  an  important  link 
is  missing  from  the  chain.  Indeed,  since  the  Cauchy  problem  for 
(3.5),  with  initial  data  of  small  total  variation,  has  global 
solutions  of  class  BV,  it  is  natural  to  expect  the  same  from  (3.1). 
This  conjecture,  however,  has  not  yet  been  verified.  The  difficulty 
lies  in  that  the  existence  theory  for  (3.5)  has  been  established  by 
means  of  Glimm's  scheme  [8]  which  hinges  on  the  explicit  construction 
of  solutions  to  the  Riemann  problem  for  systems  of  hyperbolic  con¬ 
servation  laws.  Numerical  analysts  [2]  have  applied  Glimm  type 
methods  (in  combination  with  fractional  steps  or  by  modifying  the 
Riemann  problem)  to  systems  containing  "parabolic"  terms  with 
excellent  results.  Nevertheless,  the  theoretical  justification  of 
this  approach  is  still  lacking. 

Solutions  of  (3.1)  of  class  BV  have  different  geometric 
structure  than  corresponding  solutions  of  (3.5).  For  functions 
(u,v)  of  class  BV,  first  derivatives  ut»ux,vt,vx  are  l°cally 
finite  Borel  measures.  When  such  a  pair  of  functions  is  a  solution 
to  (3.5),  the  set  of  points  of  jump  discontinuity  is  (cf.  [7))  the 
countable  union  of  C*  curves  (shocks)  with  slope 

(3.7)  s  .  »<-  Ifip)1/Z. 

The  jumps  of  u  and  v  across  shocks  are  controlled  by  the  Rankine- 
Hugoniot  conditions 


w 


(3.8) 


s[u]  +  [v]  =  0 
S[v]  -  [p(u)  ]  -  0. 
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On  the  other  hand,  if  (u,v)  e  BV  is  a  solution  to  (3.1),  it 
follows  from  (3.1)2  that  vxx  is  also  a  locally  finite  Borel 
measure.  Thus,  v  cannot  sustain  any  jump  discontinuity  and,  as 
seen  from  (3.1)^,  u  may  only  be  discontinuous  along  stationary 
lines,  s  =  0.  Through  explicit  construction  of  particular  solutions 
of  (3.1)  one  shows  that  such  singularities  may  indeed  occur.  Con¬ 
sequently,  in  contrast  to  the  complete  parabolic  damping  of 
Section  2,  viscous  damping  is  incapable  of  smoothening  rough  initial 
data. 

My  conjecture  is  that  System  (3.2)  exhibits  exactly  the  same 
behavior  as  (3.1).  A  research  program  towards  verifying  this  con¬ 
jecture  is  in  progress. 

In  System  (3.3)  damping  is  weaker  than  in  (3.1)  or  (3.2). 

Slemrod  [14]  has  shown  that  global  smooth  solutions  exist  j 

provided  that  the  initial  data  are  both  smooth  and  "small".  My  ! 

conjecture  is  that  smooth  waves  of  large  amplitude  break.  Since  j 

the  Cauchy  problem  for  (3.6)  has  a  global  solution  of  class  BV  when 
the  initial  data  have  small  total  variation  [11],  it  is  natural  j 

to  expect  similar  behavior  from  (3.3).  However,  this  conjecture  has 

not  been  established  yet.  j 

1 

Jump  discontinuities  of  BV  solutions  of  (3.6)  occur  across  ] 


forward  or  backward  shocks,  propagating  with  speeds 
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(3.9)  s-  ±{.  MftP)"2, 

or  across  contact  discontinuities  which  are  stationary,  s  *  0, 
and  across  which  u  and  0  may  jump  but  v  and  p  are  continuous. 
On  the  other  hand,  if  (u,v,0)  €  BV  is  a  solution  of  (3.3),  ®xx 
is  also  a  locally  finite  Borel  measure.  It  follows  that  6  cannot 
sustain  jump  discontinuities.  Forward  and  backward  shocks, 
propagating  at  speeds  shown  in  (3.9),  are  still  possible.  Contact 
discontinuities,  however,  are  ruled  out  since  u  alone  cannot  jump 
when  all  v,9 ,  and  p(u,0)  are  continuous. 


4.  Viscous  Damping  Induced  by  Memory  Effects 

In  this  section,  we  discuss  dissipation  mechanisms  induced  by 
viscosity  in  simple  materials  with  fading  memory  (e.g.  [16]).  As 
representative  examples  consider 


(4.1) 


0 


(4.2) 


u.  -  v  +  P  f  e"*1^  T^v  dt  *  0 
t  x  j  x 

vt  ‘  °<u>x  '  \  a^t'T^(u)xdT  "  °* 


System  (4.1)  is  a  model  of  the  equation  of  motion  of  non-linear,  one 
dimensional  viscoelastic  materials  of  the  Boltzmann  type  while  (4.2) 


is  the  equation  of  motion  of  another  model  material  (v  is  velocity  in 
both  (4.1)  and(4.2)  while  u  is  a  strain  in  (4.1)  and  an  internal 
variable,  functionof  the  history  of  velocity  gradient,  in  (4.2).  The 
functions  o  and  <f>  are  smooth  and  strictly  increasing  and  the 
kernel  a(t)  is  a  relaxation  function. 

Let  us  generally  consider  systems 

(4.3)  wt  +  f(w)x  +  |  A(t-T)g(w)xdf  =  0 

where  f  and  g  are  smooth  maps  from  Rn  to  Rn  and  A(t)  is 
a  smooth  n  x  n  matrix-valued  relaxation  kernel.  We  shall  be  look¬ 
ing  for  solutions  near  w  =  0  and  we  will  be  assuming  that  Vf(0) 
has  n  real  distinct  eigenvalues  (strict  hyperbolicity) ,  none  of 
them  zero,  in  order  to  exclude  the  possibility  of  stationary  waves  that 
would  overburden  the  memory  term. 

In  systems  of  the  above  form,  dissipation  may  be  induced 
by  the  memory  term.  In  order  to  expose  the  instantaneous 
component  of  damping,  we  define  a  new  relaxation  function 

(4.4)  B (t)  =  A(t)Vg(0)[7f(S))'1 

•V 

and  write 

(4.5)  A(t)g(w)  -  B(t)f(w)  +  A(t)h(w) 

**»  %  "»  *<» 


where 
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(4.6)  h(w)  -  g(w)  -  7g(0)[Vf(0)J_1f(w). 

Then  (4.3)  takes  the  form 

(4.7)  wt  +  f(w)x  +  J*  B(t-t)f(w)xdx  =  -  Jt^>A(t ■  x) h (w) xdt . 

Next,  we  consider  the  resolvent  kernel  K(t)  of  B(t),  i.e.,  the 
solution  to  the  linear  Volterra  equation 

(4.8)  K(t)  +  ftB(t-x)K(x)dx  =  -B (t) . 

J0 

Forming  the  convolution  of  (4.7)  with  K(t)  and  after  a  simple 
computation  we  arrive  at 

(4.9)  wt  +  f(w)x  +  K(Q)w 

=  -j  K' (t-x)wdx  +  |  E(t-x)h(w)xdx 

where 

(4.10)  E(t)  -  A(t)  +  [  K(t-x)A(x)dx. 

Equation  (4.9)  is  more  convenient  than  the  original  form,  (4.3). 

In  the  first  place,  if  B(0)  is  negative  definite,  then  the  term 
K(0)w  on  the  left-hand  side  of  (4.9)  induces  instantaneous  damping. 
Furthermore,  the  integral  terms  on  the  right-hand  side  of  (4.9)  are 
tame,  the  first  one  because  it  is  linear  in  w  and  does  not  involve 
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any  derivatives  and  the  second  because  it  is  "small",  by  virtue  of 
Vh(0)  -  0. 

The  damping  in  (4.9)  is  quite  weak.  Even  so,  in  several 
special  cases,  including  (4.1)  with  an  appropriate  relaxation 
kernel,  it  has  been  shown  [4,5,12,13,15]  that  there  exists  a  global 
smooth  solution  to  the  Cauchy  problem  provided  that  the  initial  data 
are  both  smooth  and  "small".  On  the  other  hand,  when  the  initial 
data  have  small  total  variation,  L.  Hsiao  and  the  author  [6]  have 
established,  by  means  of  a  modification  of  the  method  of  Glimm, 
the  existence  of  a  global  solution  of  class  BV.  The  strategy  in 
[6]  is  to  show  that  damping  counterbalances  the  effect  of  the 
integral  terms  on  the  right-hand  side  of  (4.9). 

Clearly,  a  lot  of  work  is  still  needed  in  order  to  complete 
the  program  of  classification  of  conservation  laws  with  dissipation. 
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